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We consider the effects of a shear on velocity fluctuations in a flow. The shear gives rise to a 
transient amplification that not only influences the amplitude of perturbations but also their time 
correlations. We show that, in the presence of white noise, time correlations of transversal velocity 
components are exponential and that correlations of the longitudinal components are exponential 
with an algebraic prefactor. Cross correlations between transversal and downstream components 
are strongly asymmetric and provide a clear indication of non-normal amplification. We suggest 
experimental tests of our predictions. 
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I. INTRODUCTION 



There has been a resurgence of theoretical interest in 
shear-driven flows, ubiquitous both in nature and in the 
laboratory [1], because the linearization of the Navier- 
Stokes equation about a laminar shear profile gives rise 
to a non-normal linear operator. The eigenstates of this 
operator are not orthogonal so, even if the laminar pro- 
file is asymptotically stable, perturbations can grow as 
transients for a while before submitting eventually to 
viscous damping [1-5]. This effect occurs generically if 
non-normal linear operators arise during a linear-stability 
analysis and it has been discussed at length in various 
contexts [6], including the hydrodynamical one we con- 
centrate on here. 

Several authors have advocated the use of peudo- 
spectra in problems involving non-normal operators [4,7]. 
A complementary approach is to characterise the be- 
haviour of non-normal systems by studying their re- 
sponse to an externally imposed noise. Farrell and loan- 
nou [8,9] showed that, in general, this results in an in- 
creased variance; they determined amplification rates 
and the perturbations that give rise to the largest am- 
plitudes for shear flows [10]. Rather little has been said 
about temporal correlations in such systems; an impor- 
tant, early exception is the study by Onuki [11] whose 
emphasis is quite different from ours as we explain be- 
low. But the temporal relation between lift-up and other 
instabilities is at the heart of the steady cycle that Wal- 
effe et al propose as the main mechanism for sustaining 
turbulent fluctuations in shear flows [12-14]. The idea is 
that downstream vortices give rise to downstream streaks 
by non-normal amplification, that the streaks go unsta- 
ble to vertical vortex formation and that finally these 
vertical vortices are folded over into downstream vortices 
by the linear shear profile. 

It is our aim here to present results on correlations in 
linear, non-normal systems driven by white noise and to 



discuss their relevance to perturbed shear flows. Specif- 
ically, we will show that temporal cross correlations be- 
tween a streak and a vortex component can provide un- 
ambiguous signs of lift-up and non-normal amplification. 

The formalism we use to study the fluctuations is sim- 
ilar to Rapid Distortion Theory or RDT [15,16]. We as- 
sume that wc have a prescribed strong shear and follow 
a perturbation in Fourier space, using Kelvin modes. In 
order to avoid technical difficulties related to advecting 
frames of reference or time-dependent wave vectors, we 
use perturbations without variation in the fiow direction. 
Nevertheless, as we will discuss later on, there is numer- 
ical evidence that this does not significantly affect the 
main conclusion that we draw about cross-correlations. 

The outline of the paper is as follows: In Section II we 
discuss a stochastic model with two degrees of freedom, 
with one vortex and one shear component, that highlights 
the non-normal coupling and the temporal correlations 
to be expected. We then turn in Section III to a full 
discussion of temporal correlations in the stochastically 
forced Navier-Stokes equation linearised about a laminar 
shear profile. Section IV contains concluding remarks, 
suggestions for experimental tests of our predictions, and 
a discussion of the relation of our work to earlier studies. 



II. A SIMPLE MODEL 

A. Vortex-Streak Coupling in Shear Flows 

Consider a linear shear profile, Uq = Sze^- Coordi- 
nates are chosen in the meteorological convention, with x 
pointing downstream, y in the spanwise direction and z 
pointing in the direction of the shear. The Navier-Stokes 
equation for the fluid velocity u linearised around this 
flow is 

atu+(u-V)Uo + (Uo-V)u=-Vp + i/Au, (1) 



1 



where p is the kinematic prcssiirc and v the kinematic vis- 
cosity of the fluid. To keep the analysis as simple as pos- 
sible we work with the Fourier modes appropriate for pe- 
riodic boundary conditions in spanwise and downstream 
directions, and free-slip boundary conditions on two par- 
allel planes in the normal direction. The analysis of the 
linear problem with Kelvin modes [17-19] shows that 
modulations in the downstream direction give rise to a 
time-dependent wave vector and faster-than-exponcntial 
damping. Farrell and loannou [10] also show that the 
most important modes for non-normal amplification do 
not have a downstream variation. Therefore, we consider 
only perturbations with wave numbers k = {0, ky, kz), 
where ky is continuous and kz — nn/d, with n an integer 
and d the distance between the bounding planes. These 
considerations help us to identify scales for the flows con- 
sidered in this section; their full significance will become 
clear in the next section. 

In order to highlight the essentials of non-normal am- 
plification, we now take a velocity field consisting of two 
modes, namely, a spanwise streak 



imax — 



(a sin ay cos Pz ' 



(2) 



(7) 



Since the maximal amplitude of the downstream compo- 
nent of the streak Us is a, the maximal modulation of 
the downstream velocity component follows to be 



i/(a2 + /32) 



Wo • 



(8) 



Note that this is proportional to the shear (the larger the 
shear the stronger the amplification) and to the inverse 
of the viscosity. This model also allows us to draw con- 
clusions about the sizes of the vortices: Assume that the 
thickness in the shear direction, i.e., the wave number /3, 
is given. Then the maximal value of Us^max is obtained 
from differentiation of Eq. (8) for a = p. li p = -K/d, 
with d the separation between the plates, the wavelength 
in the spanwise direction is about twice this width, i.e., 
A,p„„ ~ 2tt/P = 2d. 



B. Vortex-Shear Interactions with Noise 



and a downstream vortex 







Uo, = I P COS ay sin Pz 
\ —a sin ay cos Pz J 

with amplitudes s{t) and Lo{t), i.e., 

U = s{t)Us -I- Uj{t)viui ■ 



(3) 



(4) 



The two modes are sketched in Fig. 1. In the linearised 
equation the pressure disappears as the velocity fields 
are divergence- free. The; term (Uo ■ V)u drops out and 
(u- V)Uo results in a coupling between vortex and streak: 







s 

-uia'+P^) 



(5) 



The matrix on the right hand side is not symmetric be- 
cause of the coupling of both modes through the term 
(u ■ V)Uo = Su'^Bx- The dynamics that follows from the 
non-normal system (5) has exponentially decaying vor- 
tices that drive spanwise streaks: if sq and ojq denote the 
initial amplitudes, then 



s{t) = (so + Suoot)e-"^'^ +0 )*; 



(6) 



Clearly, even if there is no streak initially (i.e., sq = 
0), there will be one as time progresses as a conse- 
quence of the mixing induced by the downstream vor- 
tex. Eventually, however, both will decay. The maximal 
amplitude of the streak follows from the maximum of 
texp{—y{a^ + P^)t), which occurs at a time 



Suppose now that we add noise to the above system. 

Since the two components arc already divergence-free, 
we model the noise by adding the random perturbations 
^s{t) and ^uj{t) to each component as follows: 



-A S 
-A 



(9) 



For simplicity we assume that 



{Ut)Ut')) = DsSit-t') 

{^^{t)^u.it')) = D^5{t-t') 

{Ut)Ut')) = ^- 



(10) 



The formalism is easily expanded to more general forms. 
For notational convenience and since the quantitative de- 
pendence on wave numbers and viscosity is of minor in- 
terest in this section, we have defined A = v{a^ + P"^). If 
we start at t = — oo, the formal solution is 



u:{t)= f dt'e-^^'-^'^Ut'); 

J — oo 

s{t)= f dt'e-^^*-''ks{t') + 

J — OO 

S f dt' I dt" 

J —oo J — oo 



Ht-t')„-Ht'-t") 



un (11) 



With this choice of initial time any effects of initial con- 
ditions drop out and the effects of noise are highlighted. 
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C. Correlations 



C^,s{t,T) = {w{t)s{T)). 



(17) 



The equations for the vorticity with their exponential 

driving and a white-noise source represent an Ornstein- 
Uhlenbeck process and thus lead to exponentially decay- 
ing correlations [20,21]. Formally this follows from 

C„,^(t,r) = (^(t)a>(r)) =e-^(*+-) 

X / dt' f dT'e^^''+^'\Ut')Ur')), (12) 

J —CO J —CO 

which yields, in conjunction with Eq. (10) and the re- 
striction of the final integral up to min(t, r), 

C.,.(t,r) = |^e-^l*-l. (13) 

The correlation function for the streaks has two compo- 
nents, one resulting from the noise in the streak compo- 
nents, and one from the non-normal amplification. The 
cross terms involving the noise sources in both streak 
and vortex drop out since they are uncorrelated. The 
streak contribution is again an Ornstein-Uhlenbeck pro- 
cess with an exponentially decaying correlation function 
as in Eq. (13). The non-normal amplification gives rise 
to a contribution from the vortices, 

Cirr*'"'(^,^) = r dt' f* dt" r dr' r dr" 

J — oo J —oo J — oo J —oo 

S^e-^^'-'"+^-^"HUt")Ur"))- (14) 

Again via Eq. (10) and the constraint on the domain of 
integration to min{t',T'), the vortex contribution to the 
streak-streak correlation function becomes 

Civorte.) ^ {1 + X\t-T\) e"^!*"-! . (15) 

This shows clearly an algebraic contribution that results 
from the transient-growth characteristic of a non-normal 
system. For small time-differences, the cusp at the origin 
for the Ornstein-Uhlenbeck process becomes a rounded, 
quadratic maximum. The total correlation function 

CsAt, = + ^ (1 + ^1^1)) ' (16) 

with the time difference A = i — r, is symmetric in A 
and decays monotonically. Formally, this calculation is 
equivalent to that for a linearly damped process exposed 
to both white noise (the driving of the streaks) and ex- 
ponentially correlated noise (the cross coupling from the 
vortices) [20,21]. What is specific to our case is that 
the coloured-noise component is generated dynamically 
by the structure of the non-normal system. Cross cor- 
relations between vortex and streak can be calculated 
similarly; and they provide interesting insights into non- 
normal amplification. We define 



Then, for t > t, i.e., if the streak is probed before the 
vortex, cross correlations decay exponentially as 

a,«(i,T) = for i > r . (18) 

However, for t < t, i.e., if the streak is probed after the 
vortex, the driving of the streak by the vortex gives rise 
to the correlation function 

Co, At, r) = (1 + 2A(t - t)) e-^(^-*), (19) 

for t < T , 

which increases first, for short times, and then decreases. 
It has a maximum that lies about 21% above the value 
at zero, at A(r — t) = 1/2. The cross correlation thus 
provides a convenient and characteristic measure of the 
lift-up effect. A comparison between the three correlation 
functions discussed above is shown in Fig. 2. 

III. NOISY SHEAR FLOWS 

In order to model noise in a shear flow, we have to allow 
for spatially fluctuating velocity fields. The dynamics we 
want to study, the interaction between streaks and vor- 
tices, is mediated by the background shear flow. While 
the detailed shape of either object and also the gradient 
of the shear flow will depend on the specific boundary 
conditions, the conclusions about the correlation function 
should be less sensitive. For instance, the vortex-streak 
recycling mechanism is very similar in flows with rigid 
boundary conditions and in those with free-slip boundary 
conditions [13,14]. We, therefore, choose the theoretically 
most convenient form of free slip boundary conditions so 
that we can use a Fourier representation for the velocity 
field, 

u(x,t)=^u(k,t)e^''- (20) 

k 

and also for the noise, 

C(x,t)=^C(k,t)e*'^-. (21) 

k 

The linearized Navier-Stokes equatoin (1) becomes 

u{k, t) + Su,{k,t)e^ = -z/k2u(k,t) +$(k,t). (22) 

As in the previous section we assume that the wave vec- 
tors have no downstream component. The stochastic 
force ^ (or noise) satisfies 

UKt)}=0, (23) 
(Ci(k, t)^j (k', t')) = A{k)Pij (k)5(k + \^')S{t - t') , (24) 
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where A{k) is the ampUtude of the variance (whose pre- 
cise functional form does not matter for the moment but 
which will be discussed further at the end of this section), 
Pjj(k) = {Sij — kikj/k"^) is the transverse projector that 
enforces the incompressibility condition, and i and j de- 
note Cartesian components. We assume, for simplicity, 
that the system is infinite; the effects of boundary condi- 
tions will be included presently by suitable constraints on 
kz as in the previous section. The temporal correlation 
functions we would like to calculate (cf. Sec. II) are 

a,(t,0 = ((«i(i)%(0», (25) 

where ((• • •)) stands for an integration over all points 
in the volume under consideration and an average over 
different realizations of the noise. Given our Fourier rep- 
resentation this becomes 

Cii{t,t') = ^(ui(k,t)u,(-k,i')) ; (26) 

k 

now only an average over different realizations of the 
noise remains for every wave number k. This average 
can be calculated for each wave number individually and 
the sum can then be evaluated. Let us first write Eq. 
(22) in component form: 

u^(k, t) + SuziK t) = -uk\^{k, t) + e^(k, t); (27) 
w^(k, t) = -uk\y{k, t) + ^j,(k, i); (28) 
u^(k, t) = -vk'u,{\^, t) + ^^(k, t). (29) 

The projection onto divergence-free velocity fields and 

noise terms implies that the last two equations are cou- 
pled; they can thus be represented by a single field u^^ 
from which the y and z Cartesian components can be 
obtained as follows: 

Uy{]<i,t) ^ -u^{\s.,t)kz/k; (30) 
u^(k,t) = Ui_^{\s.,t)ky/k . (31) 

Now only two equations remain for each k: 

w^(k, t) + {Sky/k)u^{k, t) = -uk\^{k, t) + ^^(k, t); (32) 
iiuj (k, t) = -uk\uj (k, t) + ^u, (k, t) ; (33) 

and the noise terms satisfy 

(C„(k, i)6(k', t/)) = A{k)dabS{\^ + k')^(t - t') , (34) 

where the indices a and b stand for, respectively, the noise 
components in downstream (x) and perpendicular direc- 
tion (of vortex form, hence index uj). Equations (32) and 
(33) are of the same form as Eq. (9) so we can identify A 
with and S in Eq. (9) with Sky/k in Eq. (32). We 
can therefore use the correlation functions (13), (16) and 
(18, 20) calculated previously along with the expansion 
of the velocity field (25) and (26) to obtain the corre- 
lation functions between spanwise and normal-velocity 
component (subscripts 2 and 3, respectively). 



C2,2 = fcl^-(k. i)wa,(-k, t')) , (35) 

C3,3 = -^MK t)uU-K t')) , (36) 
C2,3 = Y,^^{uU^,t)u^{-k,t')), (37) 

k 

in terms of 

{u^{Kt)u^{-Kt')) = ^A{k)e-"''\'-''\ . (38) 

Thus these correlation functions decay exponentially; 
and unless there is an asymmetry in the weighting of wave 
numbers through the amplitude A{k), the cross correla- 
tion C2,3 vanishes on account of the antisymmetry under 
reflection of each wave number individually. The correla- 
tion function for the downstream component (subscript 
1) becomes 

cm(*.*')=i:^(i+IS(i+-*^i^i))«-""'^'. 

(39) 

where A = t ~ t' is the difference in times. The cross 
correlations between vortex (as measured by the normal, 
i = 3, component of the velocity) and streak (measured 
by the downstream, j = 1, component) becomes 

k 

for t>t' , 

and 

^3,i(*. t') = E + 2.fc^(*' - t))e-''''-'\ (41) 

k 

for t<t'. 

Interchanging 1 and 3 shows the symmetry Ci 3{t,t') = 

C3,l(<',t). 

The spatial dependence or, equivalently, the wave- 
vector dependence of the correlation functions is given by 
the summands in Eqs. (35)-(37). These correlation func- 
tions have high powers of the absolute value of k in the 
denominator, so all expressions are divergent for fc — > 0, 
and thus dominated by the smallest wave numbers. In 
the absence of shear, this divergence is much milder if 
we use the form of A{k) appropriate for thermal noise 
[11]. However, more singular forms have been used for 
A{k) in an effort to obtain the Kolmogorov k'^f^ energy 
spectrum at the level of a one-loop renormalization group 
[22]. Of course the nonlinear term plays a crucial role in 
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such calculations so it is not clear whether it is appro- 
priate to use such singular forms for A(k) in our study 
which does not include nonlinear ities. 

So, unless the form of A{k) favors a different selection, 
the correlations will be dominated by the smallest wave 
numbers. Given the bounding surfaces in the 2;-direction, 
separated by a distance d, we have a normal wave number 
kz = mr/d, with n a non- vanishing integer, and a con- 
tinuous spanwise wave number ky. The smallest value of 
n is n = 1 and the associated optimal ky is nn/d. Rigid 
boundary conditions may change these scales somewhat. 



IV. CONCLUDING REMARKS 

We have discussed correlations in noisy shear flow, 
both at the level of a simple model with two degrees 
of freedom and in a model with spatially varying ve- 
locity fields. The main result is that the non-normal 
coupling between downstream vortices and streaks is re- 
flected in the temporal cross correlation function between 
the downstream and normal velocity components. 

While there have been many studies of spatial corre- 
lations (Townsend [23] summarizes a large body of lit- 
erature), few have addressed temporal correlations and 
hardly any the cross correlations we described here. For 
instance, Farrell ct al. [8] study only the variance main- 
tained by forcing a noisy shear flow. Onuki [11] concen- 
trates on the case in which A{k) has the form appropriate 
for thermal noise and then calculates the renormalisation 
of the viscosity because of the nonlinear term. Bassam 
et al [24] concentrates on energy amplification but not on 
the correlation functions we consider. Blackwelder and 
Kovasznay [25] measured the cross correlations but their 
diagrams do not provide sufficient detail at short times. 

It is, therefore, only in recent numerical work on La- 
grangian models [15] that the Lagrangian version of this 
cross correlation has been considered [26]. Analysis of La- 
grangian data shows indeed the asymmetry in the cross 
correlation due to the non-normal effects. Eulerian cross 
correlations have be studied in a numerical simulation of 
turbulent shear flow [27] and they also show the cross 
correlations and their asymmetry in time [28]. 

These numerical simulations also indicate that the re- 
lation of our linear model with independent, additive and 
white in time random perturbations to natural fluctua- 
tions in a turbulent flow is not as far fetched as it may 
seem. Certain time scales and quantitative details will 
be affected by the nonlinearities neglected in our calcula- 
tion, but if there is an overall shear the main mechanism 
of non- normal growth will always be present, and that is 
at the heart of the cross correlations we study here. 

Many of the features of the cross correlations we dis- 
cuss should be accessbile experimentally. The exponen- 
tial correlation in the vortices should be reflected in 
the correlations of transversal, i.e., spanwise and nor- 
mal, velocity components. The streaks are a property of 



the downstream velocity component and the long cor- 
relations should be visible there. Finally, the streak- 
vortex correlations should be noticeable in longitudinal- 
transversal correlations. We are not aware of experimen- 
tal measurements of such correlations. The most direct 
test would involve fluctuations in laminar shear flows. 
Two experimental set iipa come to mind: 

1. One could exploit the inevitable free-stream tur- 
bulence in a wind tunnel as a noise source. A 
shear gradient can then be imposed by a flat plate, 
above which a linear shear gradient forms. Crossed 
hotwire anemometers can be used to characterise 
the statistics of the upstream noise and the veloc- 
ity field above the plate. The free decay of the 
noise between the point of characterisation and the 
point of measurement should have a negligible ef- 
fect. Reynolds numbers should be low enough so 
that no turbulence transitions in the shear layer 
are induced. 

2. There are several experiments or proposals for driv- 
ing and influencing flows in electrically conduct- 
ing fluids by electrical and magnetic fields [29-31]. 
Random perturbations of well defined strength 
could be induced by random fluctuations in the cur- 
rents. Ideally, flow geometries that are closed, such 
as a Taylor-Couettc apparatus, or spatially local- 
ized, as in a plane Couette flow [32,33] should be 
used. Such a setup should allow for a good control 
of the noise and should avoid complications arising 
from the advection of perturbations. 

Finally, we would like mention that there is growing in- 
terest in turbulent transport in magnetohydrodynamics 
in the presence of a backgroimd shear [34]. The tech- 
niques developed here can be applied fruitfully to this 
case as well, as we will show elsewhere [35]. 
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FIG. 1. Sketch of a vortex and a streak in a shear flow, 
looking downstream. The velocity components of the vortex 
lie in the plane and are indicated by the arrows. The streak 
has only one velocity component perpendicular to the plane; 
its values are indicated by a grey scale. The extrema of the 
streak (light and dark regions) are located near the maximal 
normal velocities. All scales of length and velocity are arbi- 
trary. 




FIG. 2. Correlation functions (see text) for vortex-streak 
coupling with noise of equal intensity in each component. The 
correlation functions are normalised by their values at zero, A 
is the time difference, and A, which has dimensions of inverse 
time, is deflned in the text. 
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